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Abstract. We introduce a g- weighted version of the Robinson-Schensted 
(column insertion) algorithm which is closely connected to g-Whittaker 
functions (or Macdonald polynomials with t = 0) and reduces to the 
usual Robinson-Schensted algorithm when q = 0. The g-insertion algo- 
rithm is 'randomised', or 'quantum', in the sense that when inserting a 
positive integer into a tableau, the output is a distribution of weights on 
a particular set of tableaux which includes the output which would have 
been obtained via the usual column insertion algorithm. There is also 
a notion of recording tableau in this setting. We show that the distri- 
bution of weights of the pair of tableaux obtained when one applies the 
g-insertion algorithm to a random word or permutation takes a particu- 
larly simple form and is closely related to g-Whittaker functions. In the 
case < g < 1, the g-insertion algorithm applied to a random word also 
provides a new framework for solving the g-TASEP interacting particle 
system introduced (in the language of g-bosons) by Sasamoto and Wa- 
dati SW98 and yields formulas which are equivalent to some of those 
recently obtained by Borodin and Corwin BC11 via a stochastic evo- 
lution on discrete Gelfand-Tsetlin patterns (or semistandard tableaux) 
which is coupled to the g-TASEP process. We show that the sequence 
of P-tableaux obtained when one applies the g-insertion algorithm to a 
random word defines another, quite different, evolution on semistandard 
tableaux which is also coupled to the g-TASEP process. 

1. Introduction 

We introduce a ^-weighted version of the Robinson-Schensted (column 
insertion) algorithm which is closely connected to g-Whittaker functions (or 
Macdonald polynomials with t = 0) and reduces to the usual Robinson- 
Schensted algorithm when q = 0. The insertion algorithm is 'randomised', 
or 'quantum', in the sense that when inserting a positive integer into a 
tableau, the output is a distribution of weights on a particular set of tableau 
which includes the output which would have been obtained via the usual 
column insertion algorithm. As such, it is similar to the quantum insertion 
algorithm introduced by Date, Jimbo and Miwa [DJM90J (see also |B12j) 
but with different weights. There is also a notion of recording tableau in this 
setting. We show that the distribution of weights of the pair of tableaux 
obtained when one applies the insertion algorithm to a random word or 
permutation takes a particularly simple form and is closely related to q- 
Whittaker functions. These are functions defined on integer partitions which 
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are eigenfunctions the relativistic Toda chain [Rui90, Rui99, Eti99, GLO10J 
and simply related to Macdonald polynomials (as a function of the index) 
with the parameter t = [GLOllJ. When q = 0, they are given by Schur 
polynomials. Our main result provides a starting point for developing a 
new combinatorial framework for g-Whittaker functions and related objects, 
such as Demazure and Kirillov-Reshetikhin crystals. It will be interesting 
to understand the relation to recent developments in this area, see [HHL05, 
L09, RY11, BBL11, LL12, ST12, BF12, LS13] and references therein. 

In the case < q < 1, the ^-insertion algorithm applied to a random 
word also provides a new framework for solving the g-TASEP interacting 
particle system introduced (in the language of (/-bosons) by Sasamoto and 
Wadati [SW98] and yields formulas which are equivalent to some of those 
recently obtained by Borodin and Corwin [BCllJ via a stochastic evolu- 
tion on discrete Gelfand-Tsetlin patterns — or, equivalently, semistandard 
tableaux — which is coupled to the g-TASEP process. We show that the se- 
quence of P-tableaux obtained when one applies the g-insertion algorithm to 
a random word defines another, quite different, evolution on semistandard 
tableaux which is also coupled to the g-TASEP process (after Poissonisa- 
tion). The g-TASEP process is a particular case of the totally asymmetric 
zero-range process [BKS12J. See also [BCS12J for related recent work. 

When q — > 1, the g-Whittaker functions converge with appropriate rescal- 
ing to 3 [/-Whit taker functions [GL012]. The main result of the present pa- 
per can be regarded as a natural (yet non-obvious) discretisation, in time and 
space, of the main result of the paper [012aJ, which relates a continuous-time 
version of the geometric Robinson-Schensted-Knuth (RSK) correspondence 
introduced by A.N. Kirillov [KirOlaJ, with Brownian motion as input, to the 
open quantum Toda chain with I particles. A discrete time version of that 
result has been developed in the papers [COSZ11, OSZ12J, which is formu- 
lated directly in the context of Kirillov's geometric RSK correspondence. 
The present work differs significantly from [012a, COSZ11, OSZ12J in that 
the analogue of the RSK mapping we consider here is (necessarily) ran- 
domised. In the above scaling limit, the g-insertion algorithm we introduce 
in this paper should converge in an appropriate sense to the continuous-time 
version of the geometric RSK mapping considered in [012aj, which is deter- 
ministic, and the main result of this paper should rescale to the main result 
of [012a]. This can be seen by comparing with the corresponding scaling 
limits considered in [BC11, GL012]. 

The outline of the paper is as follows. In the next section we give some 
background on the Robinson-Schensted algorithm. In Section 3, we describe 
the g-weighted version of this algorithm. The main result is presented in 
Section 4. In Section 5, we consider the g-insertion algorithm with < q < 
1 applied to a random word and explain the connection to the g-TASEP 
interacting particle system. In Section 6 we consider the algorithm applied 
to a random permutation. The proofs are given in Section 7. 
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2. The Robinson-Schensted algorithm 

The Robinson-Schensted algorithm is a combinatorial algorithm which 
plays a fundamental role in the theory of Young tableaux [Rob38, Sch61, 
Ful97, SagOO, StaOlJ. There are two versions, which are in some sense dual 
to each other, defined via insertion (or 'bumping') algorithms known as 
row insertion and column insertion. The column insertion algorithm is also 
sometimes referred to as the dual RSK algorithm, because it has a natural 
extension to zero-one matrices which was introduced by Knuth [Knu70 j . It is 
the column insertion version which we consider and generalise in this paper. 

A tableau P is a Young diagram with positive integer entries which are 
weakly increasing in each row and strictly increasing in each column. The 
corresponding diagram represents an integer partition which is referred to 
as the shape of the tableau P and denoted by shP. For example, 

112 3 
2 3 3 
3 

is a tableau with shape (4, 3, 1). To insert a positive integer k into a tableau 
P, we begin by trying to place that integer at the bottom of the first column 
of P. If the result is a tableau, we are done. Otherwise, it bumps the smallest 
entry in that column which is larger than or equal to k. Now proceed by 
inserting the bumped entry into the second column according to the same 
rule, and so on, until we have placed a bumped entry at the bottom of 
column (or on its own in a new column). For example, if we insert the 
number 2 into the tableau shown above, the outcome is 

112 3 3 

2 2 3 

3 

In this example, the 2 in the first column is bumped into the second, the 2 
in the second is bumped into the third, the 3 in the third column is bumped 
into the fourth, and the 3 in the fourth is bumped into a new fifth column 
on its own. Actually, it will be helpful for later reference to summarise this 
sequence of events in the following way: in this example, a 2 is inserted into 
the second row, and a 3 is bumped from the second row and inserted into 
the first row. 

Now, applying this insertion algorithm recursively to a word w = w\ . . . w n 6 
[l] n , starting with an empty tableau and successively inserting the numbers 
w\, u>2, • • • , w n , gives rise to a sequence of tableau P(l), P(2), . . . , P(n) = P. 
Note that it is not possible in general to recover the word w from the tableau 
P. This motivates the notion of a recording tableau, which we denote by Q. 
The tableau Q has size n and is standard, that is, it contains each of the 
numbers 1,2, ... ,n exactly once. If we denote by Q l the sub-tableau of Q 
consisting only of those entries which are not greated than i, then Q is de- 
fined by the requirement that shQ 1 = shP(i) for 1 < i < n. For example, if 
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w = 1143232 then 

1 1 3 4 1 2 5 7 

P = 2 2 g=34 
3 6 

The mapping uu \— > (P, Q) defines a bijection from the set of words [l] n to the 
set of pairs (P, Q) G 71 x S n such that shP = shQ, where 71 denotes the set 
of tableaux with entries from [I] and S n denotes the set of standard tableaux 
of size n. It is the column insertion version of the Robinson- Schensted cor- 
respondence. 

As a warm up for next section, we note that the above column insertion 
algorithm can also be described in terms of lattice paths, as follows. Suppose 
we are inserting a number k with 1 < k < I into a tableau P € 71, with 
resulting tableau P. For 1 < i < Z, set X 1 = shP 1 , and A 1 = shP\ Let 
(e,, 1 < i < /) denote the standard basis in Z'. Then A* = A* + where 
fc = > jk > ■ ■ ■ > ji > 1 is a weakly decreasing sequence defined by 

ji = max{{2 < m < j,^ : X l m \ - > 0} U {1}}, i = k, k + 1, . . . , I. 

The sequence k > jk > ife+i > •••J; > 1 determines a down/right lattice 
path in 1? from (k, k) to (Z + l,ji) by specifying the y-coordinates at which 
the path moves to the right. From the definition, this path takes a horizontal 
step to the right — > (i + l,j) whenever X l ~_\ > A* or j = 1, otherwise 
it takes a step down — > — 1). We will refer to this lattice path as 
the insertion path. The interpretation is as follows. A horizontal portion 
of the path starting at represents inserting an i into the jth row. A 
vertical portion starting at (i,j) and ending at (i,j — r) indicates that an i 
is bumped from the jth row to the (J — r)th row. For example, the insertion 
path corresponding to the previous example of inserting a 2 into the tableau 

112 3 
2 3 3 
3 

with I = 3 is illustrated in Figure 1. 

















— > 



12 3 4 



Figure 1 . An insertion path 
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3. The (/-weighted version 

In this paper, we consider the following generalisation of the column in- 
sertion algorithm. It is defined by a collection of kernels /^(P, P) which 
depend on a complex parameter q. We assume throughout that q is not 
a root of unity. If < q < 1, we interpret the quantity Ik(P, P) as the 
probability that, when we insert k into the tableau P, the output is P. 
Recall that the type of a tableaux P, which we denote tyP, is the compo- 
sition n = (/Ui,/i2, • • •) where ^ is the number of i's in P. The set of P 
for which Jfc(P, P) 7^ has the following properties. The type of P is given 
by tyP = tyP + e^. The shape of P satisfies shP = shP + ej for some 
1 < i < Moreover, if we set A 2 = shP* and A* = shP 1 , then there is a 
weakly decreasing sequence k = j^-i > jk > Jfe+i > '''iz > 1 such that 
A i = A* for 1 < i < k and A* = A* + e k for < i < I. The kernel I k (P, P) is 
defined to be zero if there is no such sequence; if there is such a sequence, 
it is given as follows. Define 



f (i,j) = 1 - q^-i-^, f^ij) = * for j > 1; 

l - gVi~ A j 

/o(i,l)=/i(t,l) = l. 



and set 



Then 




if j = ji-i and i / fc; 
otherwise. 



ji-i 



a) / fc (p,p)=n /(i.ji) n (i-/(i,i)) 

i=fc y i=ii+i 
It follows easily from the definition that 

]T>(P,P) = i. 
p 

If < q < 1, then ifc(P, P) > 0. In this case, for each k and P, Ifc(P, ■) 
defines a probability distribution on Ti and we interpret (P, P) as the 
probability that, when we insert k into the tableau P, the output is P. 

The formula (1) can be interpreted in terms of insertion paths, as follows. 
The sequence k > jk > jk+i > • • • ji > 1 determines a down/right lattice 
path in Z 2 from (k,k) to the vertical boundary {(/ + 1 < j < k} by 

specifying the y-coordinates at which the path moves to the right. The edge 
weights are f(i,j) on the horizontal edge — >• (i + 1, j) and 1 — f(i,j) 
on the vertical edge (i, j) — >■ (i, j — 1), and taking a product of these weights 
along the path gives the weight i&(P, P) for the corresponding output P. 
We interpret this path as the insertion path associated with (/-inserting the 
number k into P with resulting tableau P. As before, a horizontal portion 



6 



NEIL O'CONNELL AND YUCHEN PEI 



of the path starting at represents inserting an i into the jth row. 

A vertical portion starting at (i, j) and ending at — r) indicates that 
an i is bumped from the jth row to the (j — r)th row. When q = 0, 
there is only one output tableau P with non-zero weight, namely the output 
of the usual column insertion algorithm. Moreover, if we denote by ujq 
the insertion path corresponding to this tableau and by S(k, P) the set of 
insertion paths corresponding to the support of Ik(P, •) for nonzero q, then 
ojo £ S(k,P) and it is the 'highest' path in S(k,P) in the sense that the 
sequence k > jk > jk+i > • • • ji > 1 is maximal (in the second example 
below, it is the path shown on the top left of Figure 2). 

Let us compute the kernel Ik(P,P) for some concrete examples. 

Example 1. Suppose 1 = 2. If we are inserting a 1 into P € T2 there is 
only one possible outcome P with I\(P,P) ^ 0, namely the one obtained 
by the usual column insertion algorithm: the 1 is inserted into the first row, 
pushing the existing first row over by one. The weighted insertion path in 
this case is very simple: 





1 


1 






— > 



For example, if 



P = 



112 2 
2 



then, setting 



P^ = 



1112 2 
2 



we have 



h(P,P) 



1 if P = Pi 
otherwise. 



On the other hand, if we are inserting a 2 there are two possibilities: 



(1) The 2 is inserted into the second row, pushing the existing second 
row over by one: this outcome has weight 1 — g A i _A 2. 

(2) The 2 is inserted into the first row, pushing the existing 2's over by 
one: this outcome has weight g A i~ A 2. 
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Note that these weights sum to one, as is always the case. The corresponding 
insertion paths, with edge weights indicated, are: 



1 - q A i~ A 2 

— h 2 



1 2 

(1) 



1 2 
(2) 



The quantity A \ — A| is the difference between the number of l's in the first 
row and the number of 2's in the second row, see Figure 2. 



1 




1 


2 




2 


2 




2 


\\ — A2 

< X 





Figure 2. The quantity A} — A| in the exponent in Example 1 
For example, inserting a 2 into 



gives 



where 



and 



112 2 

2 



h{P,P)= { 



P> 







if P = P 2 
if P = Ps 
otherwise. 



112 2 
2 2 



~ 112 2 2 
P 3 = 2 

Example 2. Suppose 1 = 3. If we are inserting a 1 into P € T3 there is 
only one possible outcome P with I\(P,P) ^ 0, namely the one obtained 
by the usual column insertion algorithm: the 1 is inserted into the first 
row, pushing the existing first row over by one. The corresponding weighted 
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insertion path is: 













1 


1 


1 








— > 



12 3 4 

If we are inserting a 2, there are three possible outcomes: 

(1) The 2 is inserted into the second row, pushing the existing second 
row over by one: this outcome has weight 

/-, a 1 a 2 n i - q 1 2 

1 — g A l _A 2 

(2) The 2 is inserted into the second row, bumping a 3 into the first row: 
this outcome has weight 



(1 



1 



1 — q^i~^2 



(3) The 2 is inserted into the first row, pushing existing 2's and 3's in 
first row over by one: this outcome has weight g A i~ A 2. 
The corresponding insertion paths, with edge weights indicated, are: 



i 1 



2 

(1) 
3 

2 

1 







1 


-1*- 


X 2 




1 

i-q 


A 1 A 2 

\'i i 2 
A l A 2 













2 

(2) 













A 1 X 2 
A 1 A 2 


1 


1 








> 



2 3 4 
(3) 
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If we are inserting a 3, there are also three possible outcomes: the 3 is 
placed in the third, second or first row with respective weights 1 — g A 2 _A l ) 
qr A 2 -A i(l — and g A 2 _A lq A i _A 2. The corresponding insertion paths, 

with edge weights indicated, are: 



„A ? -A ; 

-fc 3 







X 2 X 3 
A 2 — A 3 








1 

















X 2 X 3 








X 2 X 3 


h- 1 

r 








7 







1 



The quantities \\ — A|, Af — Af, etc. which appear in the above weights 
are illustrated in Figure 3, 



a; 



A 2 



Aj — A| 



A? 



Figure 3. The exponent quantities in Example 2, 
Example 3. Suppose we are inserting a 3 into 

(2) P 



1 2 2 2 3 5 

2 3 4 5 

3 4 
5 



The (four) possible output tableaux P and their weights Is(P, P) are shown 
in Figure 4, along with the corresponding weighted insertion paths. 
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P 



1 2 2 2 3 5 

2 3 3 4 5 

3 4 
5 



I 3 (P,P) = (l-q 



1 — q 3 1 — q 2 ' 



P 



1 2 2 2 3 5 5 

2 3 3 4 

3 4 ' 
5 

1 2 2 2 3 4 5 

2 3 3 5 

3 4 
5 

1 2 2 2 3 3 5 
5_ 2 3 4 5 
3 4 
5 



I 3 (P,P) = (l-q 



1-q 3 



I 3 (P,P) = (l-q 2 )[l 



1-q 2 
1-q 3 



I 3 (P,P) = q 2 . 



1 


1-g 2 


1-^ 


1-9 



















1— 1 


l-? 2 






1 - 


1-9 2 




1 








> 



6 °2 
3 

2 

1 



I— 1 


1-g 2 


1-^ 






1 - 




1 

>— — > 











6 °2 



1 








9 2 


1 


1— 1 


1 








— > 



6 



Figure 4. The four possible output tableaux P, their 
weights I 3 (P, P), and the corresponding insertion paths, with 
edge weights indicated, for k = 3 and P given by (2). 



The ^-insertion algorithm can be applied to a word w = w\ . . . w n £ 
[l] n , starting with an empty tableau and successively inserting the num- 
bers wi, W2, ■ ■ • , w n , multiplying the weights along each possible sequence 
of output tableaux P(l), . . . , P(n) = P to obtain a distribution of weights 
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<t>w(P,Q) on Ti x S n . More precisely, we define (j) w (P,Q) recursively as fol- 
lows. Set 



1 if P = k and Q = 1 
otherwise. 



For w G [l] n and (P, Q) G 71 X 5 n+ i with shP = shQ, define 

<^*(p, Q) = 1] M^, Q)4(^, P), 

where the sum is over (P, Q) € 71 x 5 n with shP = shQ. 

We conclude this section by giving a more algorithmic description of the 
^-insertion algorithm. For this it is convenient to assume < q < 1 and 
describe it using probabilistic language, although it will be clear how to 
modify this using the language of 'weights' in the general case. For refer- 
ence, we begin with an algorithmic description of the usual column insertion 
algorithm. Denote the input word by w G [l] n . 

(1) Set *<- 1 and (P,Q) = (0,0). 

(2) Set k <— u>i and j <— k. 

(3) If ^jZl = A| and j > 1 then set j <— j — 1; otherwise k displaces the 
first number s in jth row of the tableau that is larger than k (s = oo 
and k is appended at the end of the row if no such number exists) 
and set k <— s. 

(4) If k = oo then append i to Q such that P and Q have the same 
shape, set i <— i + 1 and go to step (2); otherwise go to step (3). 

The q- insertion algorithim is defined as follows. We adopt here the following 
convention: for i > 0, let 

q A o _A i = g A o~ A i = q\)-\) = q\-\+i = q\~\ = q\ = o. 

This convention is used for covering boundary conditions in general argu- 
ments. It is only used in the following description of the g-insertion algorithm 
as well as in Section 7.1. Otherwise the undefined A*- for j > i or j = are 
taken to be zero. 

(1) Set *<- 1 and (P,Q) = (0,0). 

(2) Set k <— Wi, j <— k, d <— and a e (m, n) <— fe{m, n) Ve G {0, 1}, 1 < 
n < m. 

(3) With probability 1 — ad(k,j) set j j — 1 and d <— 0; otherwise 
displaces the first number s in jth row of the tableau that is larger 
than k (s = oo and append k at the end of jth row if no such number 
exists) and set k <— s and d <— 1. 

(4) If A; = oo then append i to Q such that P and Q have the same 
shape, set i i + 1 and go to step (2); otherwise go to step (3). 

As is obvious, when q = it reduces to the usual column insertion algorithm. 
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4. Main result 

The weights (p w (P,Q) are quite complicated. The main result of this 
paper is that a remarkable simplification occurs when we average over the 
set of words. Before stating the result, we first introduce two more func- 
tions on tableaux and explain their connection to g-Whittaker functions and 
Macdonald polynomials. Denote the g-Pochhammer symbol by 



(n) q ■.= (q;q) n = (l-q)...(l-q n ), 
with the conventions (n)o = (0) q = 1, and the g-binomial coefficients by 



in 



q (m) q (n - m) q 
For P € 71 with shP* = A*, 1 < i < I, writing A = X 1 , define 



«(P) 



I\ l j=2 rii=l(^i - K+l)q 

nKm = i(^-AS) g (A^ +i -Ai) 



aka)- 1 n 

i<i<*<' 



j ~ A i+i 



where 

l-i 

A l (\)=n(\ t -\ t+1 ) 
i=i 

For Q € <S n with shQ 1 = /j, 1 , 1 < i < n, define 



p{q)= n (i-^-"5+i). 



i<i<r- m}-^ _1 =i 



The functions k and p are simply related as follows. Suppose that I > n and 
P has distinct entries %\ < ii < ■ ■ ■ < i n . Denote by P € S n the standard 
tableau obtained by replacing the entry by k, for each k = 1, . . . , n. Then 



(3) 



K(P) 



P{P) 



(1-?)»A,(A)' 



Indeed, using the simple identities, 



(4) 



o 




J <? 



a 
1 



1 -q a 

T^7' 
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we have 



l<i<j<l 
\ i —X i ~ 

3 3 



n 

1 ±<i<j<l 



X) - Xf\ 



A,(A) 



-i 



n 



n 



l-q 



l<i<j<l 

v.-x 1 .- 

3 3 



l<i=j<l 

3 3 



P(P) 



(l-g)»A,(A)' 



The functions k and p are closely related to g-Whittaker functions [Rui90, 
Eti99, GLO10, GL012J. Denote by Q l the set of partitions with at most I 
parts. The g-Whittaker function with parameter a G C ! is a function on Q l 
defined by 



(5) 



PG7J: shP=A 



In [GLOllJ it is shown that these functions are given in terms of the Mac- 
donald polynomials P x (x;q,t) as 

(6) tf«(A) = A / (A)- 1 P A (a;(?,0). 

From this it follows that 



(7) 



^ a (X) = MXy^kx^m^c 



where denote the monomial symmetric functions and 
(8) k\n(q) = A; (A) £ *(P) = £^M)i<:, 

shP=A,tyP=/x f 



where K\a(q,t) are the two- variable Kostka polynomials [M98J. We recall 
that Kv(g,0) = Kyv'(0,q) = Ky u ,(q), where K Xfl (t) = K x ^(0,t) are the 
single-variable Kostka polynomials. For an extensive survey of the various 
properties and interpretations of these polynomials, see [KirOlb]. When 
q = 0, k(P) = 1 and k Xfl {0) is equal to the Kostka number K Xfl , which is the 
number of tableaux with shape A and type \i. In this case, ^ r Q (A) is given 
by the Schur polynomial 

*a(A) = s x (a) =J2K Xfl m fl (a). 

We will also consider the following functions: 

f X (q) = £ p(0). 

QG5„: shQ=A 
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Note that f x (0) = f x , the number of standard tableaux with shape A. 
The relation between f x (q) and the Whittaker functions ^ a is given by the 
following proposition, which is a straightforward consequence of (3). Define 

A(A) = n(Ai-A m ) g , 

i=i 

where l(X) denotes the number of parts in A. 
Proposition 1. For each Ahn, 

f X (a) 

lim * (1/11 i A) = ,. J - M \^ , . 
/^oo W l ) v > n !(l _ q) n A(X) 

It follows, using 

lim s x ((l/l) l ) = f x /n\, 
that f x (q) is also given, for A h n, by 

n 

To understand this in terms of specializations, recall that the exponential 
specialization exi is the homomorphism defined on the ring of symmetric 
functions by exi(p n ) = 6 n x, where p n are the elementary power sums (see, 
for example, [StaOl, §7.8]). It follows from the above proposition (or can be 
seen directly) that 

f X (q)=nl(l-qTex 1 (P X (q,0)) . 

The q- Whittaker functions *$> a are eigenfunctions of Ruijsenaars' relativis- 
tic Toda difference operators [Rui90, Rui99, Eti99, GLO10J. In particular, 



(9) L^ a 



where L is the kernel operator defined by 

Cj(A) if fi = A + ei for some !<«</, 



l(a,m) = «J 

and 

Ci(A) 



otherwise, 



1 - q x -^+i+ 1 for 1 < i < I, 
1 for i = I. 

Our main result is the following. 

Theorem 2. Let (P,Q) € 71 x S n with shP = shQ = A. Then 

(io) Yl Mp,Q) = (^ 1 k(p)p(Q). 

we[i] n 



A q- WEIGHTED VERSION OF THE ROBINSON-SCHENSTED ALGORITHM 15 

We note the following immediate extension of this identity which is useful 
for applications. The type of a word w is the composition \x = (^i,/i2, ■ ■ ■) 
where [n is the number of i's in w. For a = (a%, . . . , o/) and \x a composition, 
write = af ...a? 1 ; for w G [l] n and P G % write a w = a ty ^ and 
a p = a tyP . Now, since (j) w (P, Q) = unless tyP = ty(iu), we can write 

(11) aW MP, Q) = {\i)~ l a P K(P)p(Q). 

we[i] n 

Summing (11) over P and Q gives 

^ ^ o tt ^(P,Q) = (A,)7 1 * (A)/ A (?). 

(^,Q)e7Ix5„:shP=shQ=A«)GW n 

Note that this implies the Cauchy-Littlewood type identity 



Ahn V i / 



Theorem 2 also yields some combinatorial formulas. 

Corollary 3. Let A, [i h n with at most I parts, and let Q be a standard 
tableau with shape A. Then 



P X (a; q,0) = ]T H Q ( 



w)a w 



and 



where 



togfl]": ty(w)=/i 



Similarly, for any fixed P G 7; mi/i shape Ahn, 

/ A (g) = £ GpH, 

i«g[Z] n 



where 



G P (w) = ^Y J MP,Q)- 



Taking P to be standard with shape Ahn, this last formula becomes 
eXl (P x (q,0)) = - V H P (a), 

n I ^ — ^ 



n! 



where the sum is over permutations and Hp(o~) indicates the function Hp 
evaluated at the word cr _1 (l) . . . o" _1 (n). 
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When q = 0, Hs(w) equals 1 if the Q-tableau obtained by applying 
the Robinson-Schensted algorithm with column insertion to w is S, and 
otherwise; similarly, Gt{w) equals 1 if the P-tableau obtained by applying 
the Robinson-Schensted algorithm with column insertion to w is T, and 
otherwise. The functions Gt and H$ thus generalise the notions of in- 
equivalence and Q-equi valence, or Knuth and dual Knuth equivalence, for 
the Robinson-Schensted algorithm with column insertion (see, for example, 
[Ful97, Chapter 2 and §A.3]). 

The key ingredient in the proof of Theorem 2 is the following intertwining 
relation. Define kernel operators K and M by 

i 

K{\, P) = a p K (P% hP=x , M(P, P) = J2 °k J *( P > P). 

k=l 

Proposition 4. The following intertwining relation holds: 
(12) KM = LK 

We remark that (12) immediately yields the eigenvalue equation (9). 



5. Stochastic evolutions 
If < q < 1 and a G R l + with ]T • cij = 1, then 

(13) aW MP, Q) = (\i)- l a P K{P)p{Q) 

w£[l] n 

defines a probability measure on 71 x S n , which can be interpreted as the dis- 
tribution of the pair of tableaux obtained when one applies the randomised 
insertion algotihm to a random word w\ . . . w n with each wi chosen inde- 
pendently at random from [I] according to the probabilities oi,...,oj. If 
we denote by C{m) the shape of the tableau obtained after inserting the 
first m entries w\ . . . w m then, given the interpretation of Q as a recording 
tableau, we conclude by summing (13) over P that the sequence of shapes 
>C(1), . . . , C{n) is distributed according to 

F(£(l) =n\.. .,C(n) = /i n ) = vP (//>(Q), 

where Q £ S n is defined by shQ* = //, i = 1, . . . , n. But this can be written 
as 

P(£(l) =//,..., C{n) = M ») = J] ^j^ Lfr^pi). 

Since n is arbitrary, we immediately conclude the following. Write \i /*■ A if 
A is obtained from \x by adding a single box. 

Theorem 5. When applying the randomised insertion algorithm to a ran- 
dom word w\W2 ■ ■ ■ with each Wi chosen independently at random from [I] 
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according to the probabilities 01, . . . , cti the sequence of tableaux V(n), n > 
obtained evolves as a Markov chain in 7z with transition probabilities 



k=l 

The sequence of shapes C(n) = shV(n) evolves as a Markov chain in O! with 
transition probabilities 

The conditional law ofV(n), given {£(1), . . . ,C(n); C{n) = A}, is 

F(V(n) = P\ £(1), . . . , £(n); C(n) = A) = 

The conditional law of tyVin), given {£(1), . . . , £(n); C{n) = A}, is 

F(tyV(n) = n\ £(1), . . . , £(n); C{n) = A) = 

T/ie distribution of C{n) is given by 

KA) := P(£(n) = A) = (AjJ-^aCAj/^g). 

The probability distribution ^ is a particular specialisation (and restric- 
tion to A h n) of the Macdonald measures introduced by Forrester and 
Rains [FR05J, see also [BC11]. When q = 0, the above theorem reduces to 
the fact [O03J that, when applying the usual column insertion algorithm to 
a random word with probabilities a\, . . . , ai, the shape of the tableau evolves 
as a Markov chain with transition probabilities 

If ai > a<i > • • • > ai this Markov chain can be interpreted as a random 
walk in with transition probabilities 

r(M,A) = a x -n^\ 

conditioned never to exit the Weyl chamber {A € : Ai > • • • > A;}, which 
can be identified with Q l . This result, which relates to the representation 
theory of qIi, has been generalised to arbitrary complex semisimple Lie al- 
gebras in [BBO05, LLP12J. For earlier related work on the asymptotics 
of longest monotone subsequences in random words, see [TW01]. When 
q — > 1 the g-Whittaker functions converge with appropriate rescaling to 
0^-Whittaker functions [GL012], and the above theorem should re-scale to 
the main result of the paper [012aJ, which relates a continuous-time version 
of the geometric RSK correspondence introduced by A.N. Kirillov [KirOlaJ, 
with Brownian motion as input, to the open quantum Toda chain with I 
particles. In this scaling limit, the q- insertion algorithm should converge in 
an appropriate sense to the continuous-time version of the geometric RSK 
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mapping considered in [012a], which is deterministic. The results of [012a; 
have been generalised in [Chhl2] (see also [BBO09J) to arbitrary complex 
semisimple Lie algebras. It is natural to expect the results of the present 
paper to admit a similar generalisation. 

Example 4. The rank-1 case (I = 2) of Theorem 5 is discussed in [012bJ. 
Setting C l {n) = shP*(n), the evolution on tableaux in this case is driven by 
the process Y{n) = C\(n) — C\{n), n > 0, which (setting p = a±) is a birth 
and death process as illustrated in Figure 5, 



(l-p)q (1- P )q 2 (1- P )q y - 1 (1 - P )q v (1 - p)q y 




(1 - p)(l - q)(l -p)(l- q 2 ) (1 - P ){1 - q\\ -p)(l - q v+1 ) 



Figure 5. The birth-and-death process Y 

Example 5. When I = 3 the algorithm is more complicated than in the 
I = 2 case because the push-or-bump probability /i(3, 2) appears. In this 
case the algorithm with random input is described as follows (cf. Example 
2). In the following, w.p. means "with probability". 

• w.p. ai, insert 1 to row 1, pushing 2's and 3's in row 1 

• w.p. ci2, insert 2 

— w.p. 1 — g A i~ A 2, the 2 is inserted to row 2 and the displaced 3 
is either pushed or bumped 

* w.p. (1 — g A i~ A 2)/(l — g A ?^ A l ) the displaced 3 is pushed 
in row 2 

* w.p. 1 - (1 - g A ?~ A 2)/(l - q x l- x l) the displaced 3 is 
bumped to row 1 

— w.p. q A i~ A 2, the 2 is inserted to row 1 and it pushes 3's in row 
1 

• w.p. as, insert 3 

— w.p. 1 — q A 2 _A l ) the 3 is inserted to row 3 

— w.p. g A 2~ A 3(l — the 3 is inserted to row 2 

— w.p. g A 2~ A 3g A i _A i the 3 is inserted to row 1 

The g-insertion algorithm applied to a random word is closely related to 
the g-TASEP interacting particle system. This is a variation of the totally 
asymmetric simple exclusion process (TASEP) which was introduced (in 
the language of g-bosons) and shown to be integrable by Sasamoto and 
Wadati [SW98J, and recently related to g-Whittaker functions by Borodin 
and Corwin [BC11]. The process is defined as follows. There are I particles 
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on the integer lattice, and we denote their positions by x\ > X2 > • • • > x\. 
Let 01, ci2, ■ ■ ■ , ai £ R+. Without loss of generality we can assume ^2% (H = 1. 
The particles jump independently to the right by 1 with respective rates 

Joi, ifi = l; 

yaiil — q^- 1 Xi ), otherwise. 

Note that when Xi + 1 = the rate vanishes, thus enforcing the ex- 
clusion rule. Now consider the tableau- valued Markov chain V(n), n > 0, 
defined as above by applying the randomised insertion algorithm applied to 
a random word with probabilities eti, . . . , aj. Setting C l (n) = shV l (n), we 
see that the process Xi(n), . . . , Xi(n), n > defined by Xi(n) = C\{n) — i+1 
evolves as a Markov chain with state space {x € l) : x\ > X2 > ■ ■ ■ > x{\ 
and transition probabilities 

tt(x, x + ei) = ri, i = 1, . . . , I n(x, x) = 1 - y] r i} 

i 

where are defined as above. In other words, it is a de-Poissonisation of 
the g-TASEP process. Denote the g-TASEP process by X(t), t > 0, started 
with step initial condition Xj(0) = 1 — i = 1, • • • , I; by Theorem 5, the law 
of the position of the last particle at time t is given by 

F(X l (t) = m-l + l)=J2^ E ( A »)« 1 *a(A)/ A (g) 

(14) =e"* E ^W 1 *.^)/^?)- 

Aen ! ,Ai=m 

In [BCllJ, a continuous-time Markov chain on the set of tableaux Ti (actually 
discrete Gelfand-Tsetlin patterns, but this is equivalent) was introduced. It 
has the same fixed time marginals as the Poissonisation of the process V(n), 
although the dynamics are quite different. It is also coupled in exactly the 
same way to the g-TASEP process and in the paper [BCllJ an equivalent 
expression to (14) is obtained via this coupling for the law of Xi(t). See 
also [BCS12J for related recent work. 



6. Permutations 
If I = n and P E S n with shP = A, then (3) becomes 

P{P) 



(l-^A n (A)' 



Using this, and the fact that <j> w (P,Q) = unless tyP = tyu;, we immedi- 
ately deduce from Theorem 2 the following corollary. 
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Corollary 6. For P,QG S n with shP = shQ = A, we have 



(15) 



(p,q(i) ■= E MP,Q) 



p{p)p{Q) 



(l-g)"A(A)- 



Summing over P and Q gives 



Cp,q{q) 



P,QeS n :shP=shQ=\ 



(l-^A(A)- 



We note that Ylxhn^^(l) = n - When < q < 1, the probability mea- 
sure on integer partitions defined by p q (\) = 6\{q)/n\ gives the law of the 
shape of the tableaux obtained when one applied the randomised inser- 
tion algorithm to a random permutation. It would be interesting to un- 
derstand the analogue in this setting of the longest increasing subsequence 
problem [AD99, BDJ99, OkoOl]. 

For any standard tableau P with entries in [n] and shape A. Its weight 
p(P) is a product of n polynomials of the form of (1 — q k ) and hence p(P) is 
divisible by (1 — q) n . On the other hand, considering the ith. and i + 1th row 
in P, each time j a box is added in ith row, a factor (1 — q ) - where d is the 
difference between length of the corresponding two rows at time j - appears 
in p{P). For this difference d to reach the value of Aj — Aj+i eventually 
(which it evidently does) all the factors (1 — g), (1 — q 2 ), . . . , (1 — g Ai_Ai + 1 ) 
must appear at least once. It follows that p(P) is also divisible by A (A). 
Thus, Cp,q(q) £ ^[g] f° r each pair (P,Q) and 9\(q) € Z[g] for each A. 

For any permutation a G S n , denote by (P(a),Q(a)) the pair of tableaux 
after column inserting a, and set F a (q) = Cp (a), Q (a)(1)- When n = 2, the 
polynomials F a (q) and 9\(q) are given by 



Fnz{q) = (l-q) 2 ; F l32 (q) = l-q; F 213 (q) = (1 + g)(l - q 2 ); 

F 2 3l(q) = l-q 2 ; F3l2(q) = l~q 2 ; F 321 (q) = (1 + q)(l + q + q 2 ). 

03(q) = (l + q)(l + q + q 2 ); 6 21 (q) = (1 - q)(2 + q) 2 ; 9 r s (q) = (1 - q) 2 . 



The polynomials F a {q) give an alternative interpretation of the proba- 
bility measure p, q as the distribution of the shape of the tableaux obtained 
when one applies the Robinson- Schensted column insertion algorithm to a 
permutation chosen at random according to the distribution F cr {q)/n\. 



F 12 (q) = l-q; 



F 21 {q) = l+q. 



2 (g) = i + g; 



ia (g) = l-g. 



When n = 3, we have 
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7. Proofs 

7.1. Proof of Proposition 4, To prove (12), we take advantage of the 
recursive structure of the g-Whittaker functions. Define k on Q l x 1 by 



£(A',A' 



n!=i(^!. 1 _ <M + i)g(A| — a| 1 ) g 

and set 

T = {(A', A'" 1 ) 6 x n 1 - 1 : A'" 1 -< A'}, 

where we write A -< fj, if fii+i < Xi < fj.i for each i. 

We begin by verifying the simpler intertwining relation: 

(16) KM = LK, 

where M : T x T — > M>o and K : T — >• M>o are defined as follows. 

l-i 

MdX^X^^+e^X 1 - 1 )) = o,(l - g A ^- A i)n^ r ~ A ' +1 ' 1 < ^ < ^ 

i=fc 



M((X l ,X l - l ),(X l +e k ,X l - l + e k )) 



1 — g A fc-l~ A fc 

1< fc < / - 1; 



M{(X\X l - l ),(X l + e k ,X l ~ l + e m )) 

= (1 ~ - g^-^ 1 ) (1 _ gA ^- Ai) ^ 

l-g^-i-^ 1 ^ 

l<A;<m</-l. 

^ ) (A' ! A'- 1 )) = a^U^-E^A' r ^ (A ^ A /-i )lAi= „ i 

With a slight abuse of notation we will write ^(A^A'" 1 ) as shorthand for 
K(X l , (X 1 , A'" 1 )) since the support of latter is in {A' = A'}. We'll do the 
same for kernel K. 

We will verify the recursive intertwining relation (16) directly. The left 
hand side is given by 

KM(X l ,(X l + e k , X 1 - 1 )) = K(X\ A'- X )M((A', A*" 1 ), (A 1 + e k , X 1 ' 1 )) 

+ K{X\ X 1 - 1 - e k )M((X l , X 1 - 1 - e k ), (X 1 + e k , X 1 ' 1 ))!^^ 
l-l 

+ ^(A i ,V- 1 -e m )M((A',A'- 1 -e m ),(A z + e fc ,A'- 1 ))W 2 . 

m=k+l 
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We calculate each term separately. Set K' = aiK(X l , X 1 ^ 1 ). 
KiX^X'-^MdX 1 , X 1 - 1 ), (X 1 + e k , X 1 - 1 )) = K'(l - g^-*"**) J] q^'^+K 



i=k 



K(X l ,X 1 - 1 - e k )M((X l , X 1 - 1 - e k ), (X 1 + e kt X l ~ 1 )) 

/ (l-g A *-i~ A * +1 )(l-g A <= _A fe+i) (1 - g A fc " A *+i)(l - 9 A fe-i- A fc) 
_ g k — q^k *k Ti j 1 — q fe-i fc T 



_ x i 1 - g A *= A fc+i 



1 _ gA« fc -A*,- A +l 

£ i^V- 1 - e m )M((A / , A'" 1 - e m ), (X 1 + e k , X 1 ' 1 )) 

m=k+l 



A - 1 - A m+iVl _ /.AL-A^+n 



m=fe+l (l-^ Am A -+l)(l- r n- 



\ Z — 1 \ I— 1 \| \ ^ — lii \ Tfl 

(l_ g A m A ra+1 )^ _ q X m -X m +1) ^ ^-i_ AiA -q r A5-J_A!. 

=fc+l 



1 _ (7 A m-l _A m + 1 



-(1 - q x k-i~ x k) J7 »" : " 



i— 1 m 

= ir'(l-g A ^- A i) £ (l-g^-^Ui) n ^ Zll_Ai - 

m=fc+l i=fc+l 

The left hand side of (16) is thus given by 
LHS = K'(l - g A l-i- A i) (j] (? A '" 1 - A Ui 

l-l m 1 A i-1 -Ai \ 

m=k+l i=k+l 1 — q k k J 

= K(l-qh-i k) -jzt— ■ 

The right hand side is much easier to calculate: 

LK(X l , (X 1 + e k , X 1 - 1 )) = L(X l ,X l + e k )K(X l + e k , X 1 ' 1 ) 

= K' 1 - ^- A fc+i +1 *- — j-j—, 

as required. 

We will now prove (12) by induction on I. When I = 2, since M 2 is 
the kernel for the whole tableau, the recursive intertwining relation (16) is 
equivalent to the full intertwining relation (12). Suppose the statement of 
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the proposition holds for the rank-(/ — 2) case, that is, for I — 1. From the 
definition of K and K we have 

(17) K l {\\ \ ld - 1 ) = K l ~\\ l ~ x , \ ia - 2 )k\\\ X 1 - 1 ). 

By the recursive nature of definition of 4> w , M l can be expressed in terms of 
M\ M 1 - 1 and L 1 ' 1 : 

M l {X ld , X ld ) = I A i_i = ^_iM'((A', A'" 1 ), (A 7 , A*- 1 )) 

For partitions A, /U write A ~» /U to mean that either A = fi or A /* ^u. Then 
K l M\X l ,{X\X vl - 1 )) = ^(A',A u - 1 )M i ((V,A u - 1 ),(A',A 1: '- 1 )) 

A 1:,_1 :A ,_1 ~*A i_1 

= ^ ^'(A 1 , A'-^^'-HA'- 1 , A 1: '" 2 ) (l^- 1=Ai -iM'((V, A^ 1 ), (A', A'- 1 )) 



-Ml 

A'- 1 ) 

: hi-^xi-i 1 + h- va'- iIL 



+ '^■ V^-,^ ) ^ x ' _1) ' <*'■ A ' _1)) ) 



n= Yl ( ^(a', A' _1 )M'((A J , A j_1 ), (A j , A j_1 )) 

2 K'- 1 ^- 1 ,^ 



L'-ifA'- 1 A'- 1 ) 
j k\x\ X l - l )M\{X\ X 1 - 1 ), (A 1 , A'" 1 )) 



induction 



X L'-^A'- 1 ,^- 1 ) 



V k\x\ X l - l )M l ((X l , X 1 - 1 ), (X 1 , A'" 1 )) 

assumption ^ — ' \ 
A i-l^, A i-l \ 

^ L 1 ^ 1 K 1 ^ 1 (X 1 ^ 1 , (X 1 ^ 1 , A 1: '~ 2 )) \ 
L'-^A'-SA'- 1 ) y 

= ^(V,A^ 1 )Af i ((A',A^ 1 ),(A',V- 1 ))^- 1 (A i - 1 ,A 1: '- 2 ) 

A i_i/, A i_i 

Due to the indicator, when A' -1 = X 1 ^ 1 , 

I = K\X l , A^-^^-^A^ 1 , A 1: '- 2 )M'((A', A*- 1 ), (A 1 , A'- 1 )). 
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Therefore 

K l M l (X l , (X l ,X ld ~ 1 )) 

^2 K\x\X l - 1 )K l -\X l -\X ld - 2 )M l ((X\X l - 1 ),(X l ,\ 1 ' 1 )) 

= K 1 - 1 ^ 1 - 1 , X 1:l ~ 2 )K l M\X\ (A 1 , A'" 1 )) 

( = j if 1 " 1 (A'" 1 , x ld - 2 )L\x\ x l )k\x\x 1 - 1 ) 

'= L'(A', X 1 )K\X\ X 11 ' 1 ) = LK(X l , (X 1 , A 1 ^ 1 )), 
as required. 

7.2. Proof of Theorem 2. We will prove the identity (11), from which the 
statement of the theorem follows. From the definition of (j> w , for (P,Q) E 
71 x S n such that shP = shQ = A and \x l = shQ 1 for i = 1, . . . , n, the left 
hand side of (11) can be written as 

w£[l] n 

= E E a™I m (<b,P(l))...I Wn (P(n-l),P) 

= E E (a Wl I Wl ((t),P(l)))...(a w J Wn (P(n-l),P)) 

«>eM n (p(j))^- 1 1 ,shP(i)=^ 

E E W(^(l)) ••• E a Wn I Wn (P(n-l),P) 

(P(i))£T 1 1 ,shP(i)= At * / \wne[i] 

^ M(0,P(l))...M(P(n-l),P). 

(p(i))ir 1 i ,shP(i)=^ 

On the right hand side, from the definition of p(Q) and the intertwining 
relation (12), 

a P K(P)^ = L(0, M 1 ) • • ■ L(/j, n ^ 1 ,X)K(X, P) 

= L(0, /z 1 ) . . . L(/z n " 2 , n^LK^- 1 , P) 
= L(0, m 1 ) . . . L(/z n - 2 , f^-^KMi^- 1 , P) 

E (L(0,/, 1 )...L( / ."- 2 ,^- 1 ) 

P(n-l): shP^-l)^"- 1 V 



x K(/i n - 1 ,P(n - l))M(P(n - 1),P) 
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P{n-l),P{n-2): V 
shP(n-l)=^t n - 1 ,shP(n-2)=/i n - 2 

x M(P(n - 2),P(n - l))M(P(ra - 1),P)J 
= ■■■= E U(0,^)K( / . 1 ,P(l))M(P(l),P(2))x 

(P©)™"/: shP(i)=M l ^ 

•••x M(P(n-l),P)Y 

Now, from the definition of L, K and M, for P(l) £ 71 that has only one 
entry and whose shape is fi 1 = (1), 

1(0, m 1 ) = 1; P(l)) = M(0, P(l)) = a k . 

This completes the proof. 

7.3. Proof of Proposition 1, Let A h n and note that, for Z > n, Aj(A) = 
A (A). We want to show that 



f x (q) 

"I'm /,/(/\( — 



Um ^)'( A ) = n , (1 _ g) n A(A) 



From the definition of this is equivalent to 



iim r n V k(P) = — — 
i->oc _ . v ' n\(l 

Write 



.. (l-q) n A(XY 
PeTr- shP=x y H; v ' 



E ^) = ^ + ^ 

PG7I: shP=A 

where A denotes the sum over tableaux with distinct entries and B denotes 
the remaining sum. Assume I > n. By (3), if P has distinct entries, then 

<p) - p{p) 



(1-«)»A(A)' 
Hence 

A p(Q) , f x (q) 



r n A = r 



as I — > oo. Thus it remains to show that I n B — > 0. We first show that 
k(P) is bounded for P € 71 with shP = A. To see this, observe that if P 
has entries from the set . . . , z m } where i\ < ■ ■ ■ < i m and P denotes the 
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tableau obtained from P by replacing ik by k, for each k = 1, . . . ,m, then 
K (P) = K (P). It follows that 

n(P) < max k(T) < oo. 
~ rer„ 

Now, by the usual Robinson-Schensted correspondence, the number of P € 
Ti with shP = A which don't have distinct entries is at most the number of 
words w € [l] n which don't have distinct entries, and this is given by 

N(l,n) = l n - 

Clearly, l~ n N(l,n) — > as / — > oo, so we are done. 
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